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The Modified Newtonian Dynamics (MOND) generically predicts a violation of the strong 
version of the equivalence principle. As a result the gravitational dynamics of a system 
depends on the external gravitational field in which the system is embedded. This so-called 
external field effect is shown to imply the existence of an anomalous quadrupolar correction, 
along the direction of the external galactic field, in the gravitational potential felt by planets in 
the Solar System. We compute this effect by a numerical integration of the MOND equation 
in the presence of an external field, and deduce the secular precession of the perihelion of 
planets induced by this effect. We find that the precession effect is rather large for outer 
gaseous planets, and in the case of Saturn is comparable to, and in some cases marginally 
excluded by published residuals of precession permitted by the best planetary ephemerides. 



1 The external field effect with MOND 

The Modified Newtonian Dynamics (MOND) has been proposed ^ as an alternative to the dark 
matter paradigm^. At the non-relativistic level, the best formulation of MOND is the modified 
Poisson equation^, 



where p is the density of ordinary (baryonic) matter, U is the gravitational potential, g = 'VU 
is the gravitational field and g = \g\ its ordinary Euclidean norm. The modification of the 
Poisson equation is encoded in the MOND function ^{y) of the single argument y = g/ao, where 
qq = 1.2 X 10~^''m/s^ denotes the MOND constant acceleration scale. The MOND function 
interpolates between the MOND regime corresponding to weak gravitational fields y = g/ao <^ 1, 
for which it behaves as /i(y) = y + o{y), and the Newtonian strong-field regime y ^ 1, where fi 
reduces to 1 so that we recover the usual Newtonian gravity. 

An important consequence of the non-linearity of Eq. (1) in the MOND regime, is that 
the gravitational dynamics of a system is influenced (besides the well-known tidal force) by 
the external gravitational environment in which the system is embedded. This is known as 
the external field effect (EFE), which has non-trivial implications for non-isolated gravitating 
systems. The EFE was conjectured to explain the dynamics of open star clusters in our galaxy^, 
since they do not show evidence of dark matter despite the involved weak internal accelerations 
(i.e. below oq). The EFE effect shows that the dynamics of these systems should actually be 




(1) 



Newtonian as a result of their immersion in the gravitational field of the Milky Way. The EFE 
is a rigorous prediction of the equation (1), and is best exemplified by the asymptotic behaviour 
of the solution of (1) far from a localised matter distribution (say, the Solar System), in the 
presence of a constant external gravitational field Qe (the field of the Milky Way). At large 
distances r = |x| — >■ oo we have^ 

J/ = «,.x + ^£^^ + o(l), (2) 
rVl + AeSin^^ WJ 

where M is the mass of the localised matter distribution, where 9 is the polar angle from the 
direction of the external field Qq, and where we denote fie = fJ-iVc) and Ac = ycA*c/A*c) with 
ye = ge/o-o and /Hg = d/i(ye)/dye- In the presence of the external field, the MOND internal 
potential u = U — ■ yi shows a Newtonian-like fall-off ~ at large distances but with an 
effective gravitational constant G/jic."' However, contrary to the Newtonian case, it exhibits a 
non-spherical deformation along the direction of the external field. The fact that the external 
field gfe does not disappear from the internal dynamics can be interpreted as a violation of the 
strong version of the equivalence principle. 



2 Abnormal influence of the Galaxy in the Solar System 

In two recent papers ^'^ it was shown that the imprint of the external galactic field on the Solar 
System (due to a violation of the strong equivalence principle) shows up not only asymptotically, 
but also in the inner regions of the system, where it may have implications for the motion of 
planets. This is somewhat unexpected because gravity is strong there (we have g 3> ao) and 
the dynamics should be Newtonian. However, because of the properties of the equation (1), the 
solution will be given by some non-local Poisson integral, and the dynamics in the strong-field 
region will be affected by the anomalous behaviour in the asymptotic weak-field region. 

We assume that the external Galactic field g^ is constant over the entire Solar System.* 
The motion of planets of the Solar System relatively to the Sun obeys the internal gravitational 
potential u defined by 

u = U - ge-x, (3) 

which is such that lim,.^oc u = 0. Contrary to what happens in the Newtonian case, the external 
field ge does not disappear from the gravitational field equation (1) and we want to investigate 
numerically its effect. The anomaly detected by a Newtonian physicist is the difference of 
internal potentials, 

Su = u — , (4) 

where denotes the ordinary Newtonian potential generated by the same ordinary matter 
distribution p, and thus solution of the Poisson equation AtiN = —'iirGp with the boundary 
condition that limr_>.oo un = 0. We neglect here the change in the matter distribution p when 
considering MOND theory instead of Newton's law. This is in general a good approximation 
because the gravitational field giving the hydrostatic equilibrium (and thus p) is strong and 
MOND effects are very small. Hence is given by the standard Poisson integral. 

A short calculation shows that the anomaly obeys the Poisson equation ASu = — 47rG/9pdm, 
where Ppdm is the density of "phantom dark matter" defined by 

Ppdm = • ixVU) , (5) 

"Recall that in the absence of the external field the MOND potential behaves like U ~ —\/GMao Inr, showing 
that there is no escape velocity from an isolated system*. However since no object is truly isolated the asymptotic 
behaviour of the potential is always given by (2), in the approximation where the external field is constant. 

''For the Milky Way field at the level of the Sun we have Qe ~ 1.9 x 10"^^° m/s^ which happens to be slightly 
above the MOND scale, i.e. rj = Qe/ao ^ 1.6. 



where we denote x = — 1. The phantom dark matter represents the mass density that 
a Newtonian physicist would attribute to dark matter. In the model '^'^ the phantom dark 
matter is interpreted as the density of polarisation of some dipolar dark matter medium and the 
coefficient x represents the "gravitational susceptibility" of this dark matter medium. 

The Poisson equation A(5« = — 47rG/9pdm is to be solved with the boundary condition that 
lim^-^oo = 0; hence the solution is given by the Poisson integral 

f d^x' , 
6u{x,t) = G _ Ppdm(x,^). (6) 

J \^ X I 

We emphasise that, contrary to the Newtonian (linear) case, the knowledge of the matter density 
distribution does not allow to obtain an analytic solution for the potential, and the solution has 
to be investigated numerically. We can check that the phantom dark matter behaves like 
when r — >■ oo, so the integral (6) is perfectly convergent. 

In the inner part of the Solar System the gravitational field is strong {g uq) thus ^ tends 
to one there, and x tends to zero. Here we adopt the extreme case where x is exactly zero in a 
neighbourhood of the origin, say for r ^ £, so that there is no phantom dark matter for r ^ e; 
for the full numerical integration later we shall still make this assumption by posing x = 
inside the Sun (in particular we shall always neglect the small MOND effect at the centre of 
the Sun where gravity is vanishingly small). If Ppdm = when r ^ e we can directly obtain 
the multipolar expansion of the anomalous term (6) about the origin by Taylor expanding the 
integrand when r = |x| ^ 0. In this way we obtain'^ 

6u = J2^^''Ql, (7) 

1=0 

where the multipole moments near the origin are given by 

Ql = gJ^^ d^xppdn^SL (^^) . (8) 

Because the integration in (8) is limited to the domain r > e and ^^(l/r) is symmetric-trace-free 

(STF) there [indeed A(l/r) = 0], we deduce that the multipole moments Ql themselves are 
STF. This can also be immediately inferred from the fact that ASu = when r ^ e, hence 
the multipole expansion (7) must be a homogeneous solution of the Laplace equation which is 
regular at the origin, and is therefore necessarily made solely of STF tensors of type x^. Hence 
we can replace in (7) by its STF projection x^. It is now clear from the non-local integral in 
(8) that the MONDian gravitational field (for r ^ tq) can influence the near-zone expansion of 
the field when r — )■ 0. An alternative expression of the multipole moments can also be proved, 
either directly or by explicit transformation of the integral (8) . We have 

QL = -m{0)Si^o + {-y{dLum, (9) 

where the Newtonian potential un and the STF derivatives of the internal potential u are to be 
evaluated at the centre of the Sun. 

The multipole expansion (7) will be valid whenever r is much less than the MOND transition 
distance for the Solar System, defined by ro = y^GM/ao with M the mass of the Sun and oq the 

°Our notation is as follows: L = ii- ■ - ii denotes a multi-index composed of I multipolar spatial indices ii , • • • ,ii 
(ranging from 1 to 3); 9z, = di-i • • - Si, is the product of I partial derivatives dt = d/dx^; = x^^ ■ ■ - a;*' is the 
product of I spatial positions a;'; similarly = n'^ • • • rf' = x^ /r'' is the product of I unit vectors n' = x*/r; the 
symmetric-tracc-frcc (STF) projection is indicated with a hat, for instance x'" = STF[a;^], and similarly for 
and Ol- In the case of summed-up (dummy) multi-indices L, we do not write the / summations from 1 to 3 over 
their indices. 
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Figure 1: Left panel: profile of Q2{r) in the Solar System, for a standard choice of function ^i(j/) [see Eq. (13a)], 
ao = 1.2 X 10~^"m.s~^ and = 1.9 x 10^^" m.s"^. The MOND transition radius is shown by a dash-dotted line 
at ro ^ 7100 AU. Right panel: zoom of the central region (r < 50 AU), where the quadrupole is almost constant. 



MOND acceleration scale. This radius corresponds to the transition region where the Newtonian 
acceleration becomes of the order of the MOND acceleration oq and therefore, MOND effects 
become dominant. We have tq 7100 AU so the results (7)-(9) hold in a large volume around 
the Sun including all the planets (recall that Neptune's orbit is at 30 AU) . 

3 Results for the induced quadrupole moment in the Solar System 

So far we have elucidated the structure of the multipole expansion of the anomaly Su near the 
origin. Next we resort to a numerical integration of the non-linear MOND equation (1) in order 
to obtain quantitative values for the multipole moments.*^ 

The Sun being assumed to be spherically symmetric, since all the multipole moments are 
induced by the presence of the external field Qe in the preferred direction e, the situation is 
axisymmetric and all the moments Ql will have their axis pointing in that direction e. Thus 
we can define some multipole coefficients Qi by posing Ql = Qie^, where denotes the STF 
part of the product of I unit vectors = • • • e^K The multipole expansion (7) reads then as 

where Pi{z) is the usual Legendre polynomial and 6 is the angle away from the Galactic di- 
rection e. Although from the previous considerations the multipole coefficients Qi should be 
approximately constant within the MOND transition radius ro, here we compute them directly 
from the numerical solution of (1) and shall obtain their dependence on r. With our definition 
the quadrupolar piece in the internal field is given by 

Su2 = ^r^ Q2ir) (^cos^ e ~ (11) 

The radial dependence of the anomaly (11) is oc and can thus be separated from a quadrupolar 

deformation due to the Sun's oblateness which decreases like oc r~'^. 

As a first result, we show in Fig. 1 the profile of the quadrupole induced by the MOND theory 
through the function Q2{'r') defined in Eq. (11). We find that this quadrupole is decreasing from 



■^Our numerical scheme is based on the very efficient integrator of elliptic equations lorene, available from 
the website http : //www . lorene . obspm . f r . 



Table 1: Numerical values of the quadrupole Q2 together with the associated dimensionless quantity 92 defined 
by Eq. (12). All values are given near the Sun. We use different choices of the function //(y) defined in Eqs. (13). 



MOND function 




i^2{y) 


m{y) 


/^cxp(y) 


MTeVeS(y) 


Q2 [s-^] 


3.8 X 10-^*^ 


2.2 X 10-^*^ 


2.1 X w-'^^ 


3.0 X 10"^'^ 


4.1 X 10"^*^ 


<12 


0.33 


0.19 


1.8 X 10-2 


0.26 


0.36 



the Sun's neighbourhood to zero, on a typical scale of 10000 astronomical units (AU). However, 
we check numerically that Q2{r) is almost constant in a large sphere surrounding the Solar 
system, as it has a relative variation lower than 10~^ within 30 AU (see the zoomed region in 
Fig. 1). We shall therefore refer to the quadrupole as a simple number, noted (52(0) or simply 
Q2, when evaluating its influence on the orbits of Solar-system planets. 

On dimensional analysis we expect that the quadrupole coefficient Q2 should scale with the 
MOND acceleration oq like 

Q2 = -q2{v), (12) 



where ro = ^JGM|aQ is the MOND transition radius and where the dimensionless coefficient 
q2 depends on the ratio r/ = ^e/^o between the external field and ao, and on the choice of the 
interpolating function /x. Our numerical results for the quadrupole are given in Table 1, for 
different coupling functions /x(j/).^ Here we consider various cases widely used in the literature: 

/Xexp(j/) = 1 - e-^ (13b) 

/XTeVeS(y)= ^^^^^^^ . (13c) 

The function [i\ has been shown to yield good fits of galactic rotation curves^; However because 
of its slow transition to the Newtonian regime it is a 'priori incompatible with Solar System 
observations. The function /X2 is generally called the "standard" choice and was used in fits^°. 
We include also the function /Xexp having an exponentially fast transition to the Newtonian 
regime. The fourth choice ^utcVcS is motivated by the TeVeS theory One should note that 
none of these functions derives from a fundamental physical principle. 

We have used several functions of type as defined in Eq. (13a). One can notice that 
the value of Q2 decreases with n, that is with a faster transition from the weak-field regime 
where )u(y) ~ y, to the strong field regime where ~ 1. We have been unable to determine 
numerically a possible limit for Q2 as n goes to infinity. 



4 Effect on the dynamics of the Solar System planets 

We investigate the consequence for the dynamics of inner planets of the Solar System of the 
presence of an abnormal quadrupole moment Q2 oriented toward the direction e of the galactic 
centre. Recall that the domain of validity of this anomaly is expected to enclose all the inner Solar 
System (for distances r < ro ~ 7100 AU), with the quadrupole coefficient being constant up to 
say 50 AU (see Fig. 1). As we have seen, the anomaly induces a perturbation on the Newtonian 
gravitational potential, namely u = un + where un = GM/r and the perturbation function 
R = 5u Ss given for the quadrupole moment by Eq. (11). 



"Note that the quadrupole coefficient Q2 is found to be always positive which corresponds to a prolate elon- 
gation along the quadrupolar axis. 



We apply the standard linear perturbation equations of celestial mechanies The unper- 
turbed Keplerian orbit of a planet around the Sun is described by six orbital elements. For 
these we adopt the semi-major axis a, the eccentricity e, the inclination / of the orbital plane, 
the mean anomaly £ defined by ^ = n{t — T) where n = 2t:/P {n is the mean motion, P is the 
orbital period and T is the instant of passage at the perihelion), the argument of the perihelion 
U! (or angular distance from ascending node to perihelion), and the longitude of the ascending 
node il. We also use the longitude of the perihelion defined hy uj = u + Q. 

The perturbation function R = Su2 is a function of the orbital elements of the unperturbed 
Keplerian ellipse, say {ca} = {a,e,I,£,uj,Q}. The perturbation equations are generated by 
the partial derivatives of the perturbation function with respect to the orbital elements, namely 
dR/dcA- We express the planet's absolute coordinates (x, y, z) (in some absolute Galilean frame) 
in terms of the orbital elements {a,e, I,£,uj,Q} by performing as usual three successive frame 
rotations with angles Q, I and uj, to arrive at the frame {u,v,w) associated with the motion, 
where (u, v) is in the orbital plane, with u in the direction of the perihelion and v oriented in 
the sense of motion at perihelion. The unperturbed coordinates of the planet in this frame are 

u = a{cosU — e), (14a) 
V = a^/l - sinU , (14b) 
■u; = 0, (14c) 

where U denotes the eccentric anomaly, related to £ by the Kepler equation £ = U — e sin U. 
The perturbation equations provide the variations of the orbital elements dc^/dt as linear com- 
binations of the partial derivatives dR/dcs of the perturbation function. We are interested only 
in secular effects, so we average in time the perturbation equations over one orbital period P. 
Denoting the time average by brackets, and transforming it to an average over the eccentric 
anomaly U, we have 

fit'^^Lr'.iUil-eeosU)'^. (15) 



\ — ; — / = — f cit — — = -— I U.U (1 — c cos u ) —— 
\dt / PJo dt 27r Jo dt 

In the following, to simplify the presentation, we shall choose the x-direction of the absolute 
Galilean frame to be the direction of the galactic centre e = Qc/gc- That is, we assume that 
the origin of the longitude of the ascending node O lies in the direction of the galactic centre. 
Furthermore, in order to make some estimate of the magnitude of the quadrupole effect, let us 
approximate the direction of the galactic centre (which is only 5.5 degrees off the plane of the 
ecliptic) as being located in the plane of the orbit; consequently we choose / = 0. In this case 
a; = a; + O is the relevant angle for the argument of the perihelion. We then find the following 
non-zero evolution equations: 



de\ 5Q2eVl - . . 

- sm(2a;} , (16a) 



di / 4ra 

V + 3e^ + 15(1 + e^) cos(2a;) , (16b) 



d^\ ^ n- — 
dt/ ~ ^ 12nL 



dL0\ 



dt / 4n 



l-F5cos(2a;) . (16c) 



We recall that uj is the azimuthal angle between the direction of the perihelion and that of the 
galactic centre (approximated to lie in the orbital plane). Of particular interest is the secular 
precession of the perihelion (dtD/dt) due to the quadrupole effect henceforth denoted by 



^'^''^^^ l + 5cos(2c:;) . (17) 



- 

4n 



Table 2: Results for the precession rates of planets A2 due to the quadrupole coefficient Q2. We use the values for 

Q2 for various MOND functions as computed in Tabic 1. Published postfit residuals of orbital precession (after 
taking into account the relativistic precession). All results are given in milli-arc-seconds per century. 



Quadrupolar precession rate A2 in mas/cy 



MOND function 


Mercury 


Venus 


Earth 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 




0.04 


0.02 


0.16 


-0.16 


-1.12 


5.39 


-10.14 


7.93 


^2(2/) 


0.02 


0.01 


0.09 


-0.09 


-0.65 


3.12 


-5.87 


4.59 


M2o(?y) 


2 X 10"''' 


10-'^ 


9 X IQ-^ 


-9 X 10-3 


-0.06 


0.3 


-0.56 


0.44 


Mexp(,!/) 


0.03 


0.02 


0.13 


-0.13 


-0.88 


4.25 


-8.01 


6.26 


MTeVes(2/) 


0.05 


0.02 


0.17 


-0.17 


-1.21 


5.81 


-10.94 


8.56 






Postfit residuals for 


A = (dw/dt) in mas/cy 








Origin 


Mercury 


Venus 


Earth 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


Pitjeva 


-3.6 ±5 


-0.4 ±0.5 


-0.2 ±0.4 


0.1 ±0.5 




-6 ±2 






Fienga et al. 


-10 ±30 


-4 ±6 


0± 0.016 


0±0.2 


142 ± 156 


-10 ±8 


± 2 • 10^ 


± 2 • lO'' 


Fienga et al. 


0.4 ±0.6 


0.2 ± 1.5 


-0.2 ±0.9 - 


0.04 ±0.15 


-41 ±42 0.15 ±0.65 







The precession is non-spherical, in the sense that it depends on the orientation of the orbit 
relative to the galactic centre through its dependence upon the perihelion's longitude Co. The 
effect scales with the inverse of the orbital frequency n = 2ti/P and therefore becomes more 
important for outer planets like Saturn than for inner planets like Mercury. This is in agreement 
with the fact that the quadrupole effect we are considering increases with the distance to the 
Sun (but of course will fall down when r becomes appreciably comparable to ro, see Fig. 1). 

Our numerical values for the quadrupole anomalous precession A2 are reported in Table 2. 
As we see the quadrupolar precession A2 is in the range of the milli-arc-second per century 
which is not negligible. In particular it becomes interestingly large for the outer gaseous planets 
of the Solar System, essentially Saturn, Uranus and Neptune. The dependence on the choice of 
the MOND function ^ is noticeable only for functions Univ) defined by (13a) with large values 
of n, where the effect decreases by a factor ~ 10 between n = 2 and n = 20. 

We then compare in Table 2 our results to the best published postfit residuals for any pos- 
sible supplementary precession of planetary orbits (after the relativistic precession has been 
duly taken into account), which have been obtained from global fits of the Solar System dynam- 
13,14,15 jj^ particular the postfit residuals obtained by the INPOP planetary ephemerides^^'^^ 
use information from the combination of very accurate tracking data of spacecrafts orbiting dif- 
ferent planets. We find that the values for A2 are smaller or much smaller than the published 
residuals except for the planets Mars and Saturn. Very interestingly, our values arc smaller or 
grossly within the range of the postfit residuals for these planets. In the case of Saturn notably, 
the constraints seem already to exclude most of our obtained values for A2 , except for MOND 
functions of the type /x„ and given by (13a) with rather large values of n. 

However let us note that the INPOP ephcmcridcs arc used to detect the presence of an 
eventual abnormal precession, not to adjust precisely the value of that precession On the 
other hand the postfit residuals are obtained by adding by hands an excess of precession for 
the planets and looking for the tolerance of the data on this excess ^^'^^ in order to really 
test the anomalous quadrupolar precession rate A2, one should consistently work in a MOND 
picture, i.e. consider also the other effects predicted by this theory, like the precession of the 
nodes, the variation of the eccentricity and the inclination, and so on — see Eqs. (16). Then 
one should perform a global fit of all these effects to the data; it is likely that in this way the 
quantitative conclusions would be different. 

Finally let us cautiously remark that MOND and more sophisticated theories such as TeVeS^^, 
which are intended to describe the weak field regime of gravity (below oq), may not be extrap- 
olated without modification to the strong field of the Solar System. For instance it has been 



argued^ that a MOND interpolating function n which performs well at fitting the rotation curves 
of galaxies is given by /Ui defined by (13a). However this function has a rather slow transition to 
the Newtonian regime, given by //i ~ 1 — when y = g/ao oo, which is already excluded 
by Solar System observations. Indeed such slow fall-off —y~^ predicts a constant supplementary 
acceleration directed toward the Sun 5g^ = oq (i.e. a "Pioneer" effect), which is ruled out be- 
cause not seen from the motion of planets. Thus it could be that the transition between MOND 
and the Newtonian regime is more complicated than what is modelled by Eq. (1). This is also 
true for the dipolar dark matter model ^'^ which may only give an effective description valid in 
the weak field limit and cannot be extrapolated as it stands to the Solar System. While looking 
at MOND-like effects in the Solar System we should keep the previous proviso in mind. The 
potential conflict we find here with the Solar System dynamics (notably with the constraints on 
the orbital precession of Saturn ^^'^^) may not necessarily invalidate those theories if they are 
not "fundamental" theories but rather "phenomenological" models only pertinent in a certain 
regime. 

In any case, further studies are to be done if one wants to obtain more stringent conclusions 
about constraints imposed by Solar-system observations onto MOND-like theories. More precise 
observations could give valuable informations about an eventual EFE due to the MOND theory 
and restrict the number of possible MOND functions that are compatible with the observations. 
More generally the influence of the Galactic field on the Solar-system dynamics through a possi- 
ble violation of the strong version of the equivalence principle (of which the EFE is a by-product 
in the case of MOND) is worth to be investigated. 
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